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ABSTRACT 

Several thin, Keplerian, sub-parsec megamaser disks have been discovered in 
the nuclei of active galaxies and used to precisely determine the mass of their host 
black holes. We show that there is an empirical hnear correlation between the disk 
radius and the black hole mass. We demonstrate that such disks are naturally 
formed by the partial capture of molecular clouds passing through the galactic 
nucleus and temporarily engulfing the central supermassive black hole. Imperfect 
cancellation of the angular momenta of the cloud material colliding after passing 
on opposite sides of the hole leads to the formation of a compact disk. The radial 
extent of the disk is determined by the efficiency of this process and the Bondi- 
Hoyle capture radius of the black hole, and naturally produces the empirical linear 
correlation of the radial extent of the maser distribution with black hole mass. 
The disk has sufficient column density to allow X-ray irradiation from the central 
source to generate physical and chemical conditions conducive to the formation 
of 22 GHz HtO masers. For initial cloud column densities < lO^^'^cm"^ the disk 
is non-self gravitating, consistent with the ordered kinematics of the edge-on 
megamaser disks; for higher cloud columns the disk would fragment and produce 
a compact stellar disk similar to that observed around Sgr A* at the galactic 
centre. 

Subject headings: accretion, accretion disks — galaxies: Seyfert — Galaxy: cen- 
ter — ISM: clouds — masers — stars: formation 
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1. Introduction 

The nuclei of some Seyfert 2 galaxies are home to powerful 22 GHz water masers located 
in a circumnuclear disk within a parsec of the central massive black hole. In almost edge-on 
systems, VLBI observations of the maser kinematics enable an accurate determination of 
the black hole mass. The archetypal system is NGC 4258, for which the black hole mass has 
been measured to 1% accuracy (e.g., Herrnstein et al. 2005; Humphreys et al. 2008). Recent 
VLBA observations of Seyfert 2 galaxies have discovered a limited number of additional 
examples with high inclination angles, allowing precise measurements of the mass of their 
host black holes as well as the physical size of the maser disks, which range from 0.1-1 pc 
(Herrnstein et al. 2008; Kuo et al. 2010). 

The physical parameters of the disks can be inferred from the conditions necessary to 
generate the masers. Collisional inversion of the 22 GHz H2O transition requires densities of 
10^-10^^ cm~^, temperatures in the range 300-1 000 K, and a sufficient column of water to 
ensure maser amplification (e.g. Neufeld & Melnick 1991). These conditions are plausibly 
produced by irradiation of a molecular disk by the central X-ray source as long as the surface 
density exceeds ~ Igcm^^ and warping exposes the disk surface to the center (Neufeld et 
al. 1994; Maloney 2002). Parsec-scale maser disks must therefore have masses > IO^Mq. 

Further insight is provided by examining the center of the Milky Way Galaxy, where 
there is strong evidence of star formation occurring in a sub-parsec scale disk 6 x 10® years 
ago (Paumard et al. 2006). Approximately 100 massive stars orbit within a few tenths of 
a parsec of the ~ 4 x IO^Mq black hole Sgr A* (Bartko et al. 2010; Lu et al. 2010; Do 
et al. 2009; see also the recent review by Genzel et al. 2010 and references cited therein). 
About 2/3 of these are locahzed in a clockwise rotating stellar disk with a wide range of 
eccentricities (Levin & Beloborodov 2003), with the remainder loosely distributed, possibly 
in a larger counter-rotating disk (Paumard et al. 2006; but see Lu et al. 2009). Stellar disks 
could be created by the tidal disruption of an inspiralling stellar cluster (Gerhard 2001; 
McMillan & Portegies Zwart 2003; Portegies Zwart et al. 2003; Kim et al. 2004; Giirkan & 
Rasio 2005), but the modelling implies that this mechanism produces a far more disordered 
stellar system than observed, and a non-existent population of massive stars shed from the 
cluster extending beyond 0.3 pc from Sgr A*; furthermore the inspiralling time scale is longer 
than the stellar ages (Paumard et al. 2006; but see Fujii et al. 2008). A more attractive 
alternative is that these disks could form "in-situ" by gravitational collapse in a disk of gas 
captured by the black hole (Levin & Beloborodov 2003; Nayakshin et al. 2007), a process 
previously considered in the context of AGN (Kolykhalov & Sunyaev 1980; Shlosman & 
Begelman 1987; Collin & Zahn 1999; Goodman 2003). 

The compactness of the stellar disks relative to molecular cloud dimensions implies that 
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the raw material for the disk is captured from a cloud as it temporarily engulfs Sgr A* while 
passing through the central parscc of the Galaxy (c.f. Sanders 1981; Bottema & Sanders 
1986); the cancellation of angular momentum of the captured cloud material that passes on 
opposite sides of the black hole naturally produces a compact, gravitationally unstable disk 
that is consistent with dimensions, mass and kinematic properties of the observed stellar disk 
(Wardle & Yusef-Zadeh 2008, hereafter Paper I; BonneU & Rice 2008; Mapelh et al. 2008; 
Alig et al. 2011). This process may also be responsible for the parsec-scale circumnuclear 
ring (Sanders 1998; Paper I). 

A link between the stellar disk at the Galactic center and AGN megamascr disks was 
suggested by Milosavljevic & Loeb (2004) who argued that the masing disks are gravitation- 
ally unstable and would eventually create stellar disks analogous to those seen around Sgr 
A*. However, the maser distribution and kinematics in the edge-on systems are very thin and 
close to keplerian, so that at least in these systems the megamaser disks are gravitationally 
stable. 

In this Letter, we show that the cloud capture scenario outlined in Paper I creates 
gravitationally stable sub-parsec disks with physical conditions conducive to the formation 
of megamasers. We begin by showing that systems with well-determined host black hole 
masses and Keplerian rotation profiles display a strong correlation between disk size and 
black hole mass (§2). In §3 we consider a simple model for the partial capture of molecular 
clouds and demonstrate that this correlation is a natural outcome of the cloud-engulfment 
scenario, and that, for reasonable parameters, initial cloud column densities ~ 10^^ cm~^ will 
create a non-self-gravitating molecular disk consistent with the observed megamaser disks. 
Our conclusions are summarised in §4. 

2. A Correlation Between Megamaser Disk Size and Black Hole Mass 

Fourteen resolved parsec-scale megamascr disks in AGN have been reported in the liter- 
ature. Of these, three disks (in NGC 1068, NGC 3079, and IC 1481) have maser distributions 
that are spatially and kinematically disordered, with rotation curves that are flatter than 
Keplerian (Greenhill et al. 1996; Kondratko et al. 2005; and Mamyoda et al. 2009, respec- 
tively) . This is suggestive of possibly clumpy material distributed within a dynamically hot 
torus with the gravitational potential dominated by the torus material or a central stellar 
cluster. Three others cases (NGC 3393, NGC 4945 and Circinus) have poorly-determined 
rotation curves so that the black hole mass and dynamical state of the disk are uncertain 
(Kondratko et al. 2008; Greenhill et al. 1997, 2003 respectively). 
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As our model for the dynamics of disk formation relies on the black hole mass dominating 
the gravitational potential on a scale several times that of the disk (see §3), we focus our 
attention on the remaining eight megamaser disks with well-determined, Keplerian profiles 
and accurate black hole masses (Herrnstein et al. 2008; Kuo et al. 2011). In each case the 
upper limits to the disk mass estimated from the errors to the fit to the rotation profile are 
a few percent of the black hole mass. Table 1 lists the the measured mass of the black hole, 
and the inner and outer radii of the disk as traced by the masers for each of these eight 
megamaser disks. For comparison we have also included the range of radii of the orbits of 
S-stars forming the compact stellar disk around Sgr A* (Lu et al. 2010). The remarkable 
similarity of size scale to the megamaser disks suggest that they may have been formed the 
same way. 

For each of these disks we plot the range of disk radii traced by megamasers against 
black hole mass (see Fig. 1). Remarkably, there is a well-defined upper envelope to the disk 
radii that scales linearly with black hole mass, given approximately by 



where the central black hole mass is M = My x 10 M©. The outer edge of the mascr 
distribution may simply trace where the surface density of the disk falls below the ~ 1 g cm^^ 
needed to shield molecules from the hard X-ray flux from the central source (Maloney 2002). 
However, there is no obvious reason why this radius should follow this trend. Our preferred 
explanation, outlined in the following section, is that the outer radii indicate a genuine 
physical truncation, and that the linear dependence on black hole mass reflects the kinematics 
of cloud capture process that forms the disk (see eq. 5). 



In Paper I we showed that partial capture of molecular clouds passing through the 
central parsec of the galaxy by the 4 x 10^ Mq black hole Sgr A* would naturally create 
a disk on the 0.1 parsec scale of the stellar disk. There we focussed on the formation of 
self-gravitating disks, which because of the short coohng time of the gas, create stars on 
a timescale comparable to a few orbital times, so that star formation would occur before 
the disk settles into a thin ordered structure; the resulting stellar orbits have a range of 
eccentricities and inclinations. 

Here we apply this scenario to the very thin and well-ordered megamaser disks, with 
two key differences from our earlier analysis. First, we explicitly follow the scaling with black 
hole mass rather than flxing it at 4 x IO^Mq appropriate for Sgr A*; second, we focus on 
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0.3 M^ pc , 




3. Disk Formation by Cloud Capture 
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Table 1. Physical Characteristics of Megamaser Disks 



Source 


BH Mass 10^ M© 


Disk Size (pc) 


NGC 1194'^ 


6.5±0.3 


0.53-1.33 


NGC 2273« 


0.75±0.04 


0.028-0.084 


UGC 3789" 


1.04±0.05 


0.084-0.30 


NGC 2960" 


1.16±0.05 


0.13-0.37 


NGC 4258'' 


3.7±0.01 


0.17-0.29 


NGC 4388'^ 


0.84±0.02 


0.24-0.29 


NGC 6264^^ 


2.91±0.04 


0.24-0.80 


NGC 6323" 


0.94±0.01 


0.13-0.30 


Sgr A*'^ 


0.43±0.02 


0.038-0.13 



^Kuo et al. 2011 

''Herrnstein et al. 2008 

'^compact stellar disk; Lu et al. 2008 
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Fig. 1. — Vertical bars indicate the range of radii of the maser distribution for pubhshed 
H2O megamaser disks with well-determined Keplerian rotation curves versus central black 
hole mass (Kuo et al. 2011; Herrnstein et al. 2005) . The extent of the compact stellar disk 
at the Galactic center (lowest-mass point; Lu et al. 2008) is plotted for comparison. The 
errors in the black hole masses and radii are a few percent. The dashed line represents an 
approximate empirical upper envelope to the outer radii given by i?max = 0.3 (M/1O''M0) pc 
(see text). 



-7- 



capture events that produce non-self-gravitating disks, which therefore settle into very thin 
Keplerian disks rather than form stars. 

Consider, then, an extended cloud with hydrogen column density A^h = -^^23 x 10^^ cm~^ 
that enters the black hole's sphere of influence with speed v — 225 V225 kms~^. On its 
passage through the region, the cloud temporarily engulfs the black hole and shocks induced 
by gravitational focussing of the material passing close to the hole lead to capture of some 
material, which circularizes. Rapid cooling of the gas leads to formation of a thin, Keplerian 
disk. 

In this scenario the disk mass is determined by the mass of gas with impact parameter 
inside the Hoyle-Lyttleton radius (Hoyle & Lyttelton 1939) 

2GM Mr 

and its size is determined by the imperfect cancellation of the angular momentum of the 
captured material that approached on opposing sides of the black hole due to the cloud 
morphology and structure. As in Paper I we characterize the uncertain capture dynamics 
using two parameters. The first of these, k, is the ratio of the captured mass to the Hoyle- 
Lyttleton estimate, allowing us to write the disk mass as 

Md^K Tibl {lAmuNn) = 1.0 x 10^ ^^^f^' . (3) 

where ttih is the mass of a hydrogen atom and we have assumed the standard helium to 
hydrogen of 0.4 by mass. We expect that k ^ 1 because bulk kinetic energy is efficiently 
lost through shocking and rapid radiative cooling during the collision (e.g. Edgar 2004, and 
references therein). 

The second parameter. A, is the average ratio of the specific angular momentum of a 
fiuid element in the resulting Keplerian disk, (GMr)^/^, to its initial angular momentum bv. 
This ignores a number of complications: the cancellation of angular momentum is sensitive 
to the distribution of inhomogcncitics in the cloud and will be highly variable; it will also 
depend on the orientation of a fiuid elements trajectory relative to the final disk plane. 
While fluid elements with the same impact parameter b but with different speciflc angular 
momentum vectors will end up on different disk radii, the convolution in mapping the initial 
impact parameter b to eventual location in the disk r also has the effect of averaging out 
these differences. This assumption allows us to relate the impact parameter of a fluid element 
in the approaching cloud, to its typical location in the disk once it has been captured and 
its motion circularized: 
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The value of A is uncertain, but recent simulations by Alig et al. (2011) suggest that A ~ 0.3 
- 0.4. 

The disk radius is determined by the initial angular momentum of fluid elements with 
impact parameter bo, for which eq Q yields 

R, = 2X\^0.31^^^ PC (5) 



-"225 



and so the mean surface density of the disk is 



Strikingly, the disk size corresponds to the upper envelope of the maser disks with the 
canonical model parameters v ~ 225kms~^ and A ~ 0.3. In other words, the kinematics 
of the model naturally produces the observed extent of the maser disks, and the empirical 
scaling with black hole mass. 

To estimate the surface density profile of the disk we note that eq Q implies that a 
thin cylindrical ring of cloud material with impact parameters in the range [6, h + ends 
up in the annulus [r, r + dr\ in the disk. Then mass conservation relates the cloud and disk 
surface densities via l-nh niXAniYj^Yi) 'ih = 2'KrTj{r)dr implying that the disk has a 1/r 
surface density profile: 

S(r) = A — ^ 1.1 "2 — 2 gem . (7) 

^0.3 "^225 '"pc 

where in the final expression Vpc is the cylindrical radius in parsecs. For this profile the disk 
mass enclosed within radius r of the BH is linearly proportional to r, and the density at the 
disk edge is 0.5 T.^. 

This simple model disk must meet three requirements in addition to its size to guarantee 
that masers may be present over the radial extent of the disk. First, the mass given by eq 
(|5| must not exceed the mass of the incoming cloud, otherwise the disk will be truncated at 
a smaller radius where the enclosed disk mass equals the cloud mass Md. 

Second, coUisional inversion of the 22 GHz H2O transition requires densities of 10''- 
10^^ cm~'^, temperatures in the range 300-1 000 K, and a sufficient column of water to ensure 
significant maser amplification (e.g. Neufeld & Melnick 1991). These conditions are produced 
by X-ray irradiation of the disk as long as warping exposes the disk surface to the central 
X-ray source and the disk column density exceeds ~ lgcm~^ (Neufeld, Maloney & Conger 
1994; Maloney 2002). For lower column densities the X-rays dissociate any molecules in the 



- 9 - 



disk, for higher column densities there is a transition from atomic gas to a water-bearing 
molecular layer layer with temperatures of ~ 400 K. Because inversion of the 22 GHz line 
can occur in this layer for a broad range of densities, the maser action is insensitive to the 
X-ray flux and disk column. Therefore we need only adopt as a requirement that the column 
density of our model disks exceed lgcm~^. Third, the disk should be gravitationally stable 
for r < R, i.e. S(r) < Csfl/nG, where Cg is the isothermal sound speed at T 400 K 
and Q is the Keplerian orbital frequency. These constraints are most stringent at Rd, and 
for our 1/r surface density profile the density at the disk edge is 0.5 S^, yielding the limits 
< Mci/(7ri?2) and Igcm-^ < ^S, < c,fi(/?rf)/7rG. 

The constraints on T.^ are mapped to the cloud column density using eq ^ with v = 
225kms~^, k = 1 and A = 0.3 and are plotted in Figure 2 as a function of black hole mass. 
To recap, the column of the incoming cloud must lie within the unshaded region to guarantee 
that the disk formed by the captured material is thick enough to allow the X-rays to provide 
the physical conditions and water column needed to produce the megamaser emission, but 
not so thick that the disk either becomes gravitationally unstable or unrealistically massive 
compared to the likely mass of the cloud. From Fig. 2, we see that these conditions are 
satisfied for cloud column densities < lO^^'^cm"^, for black hole masses approaching lO^Mo 
the disk may be truncated because the entire cloud is captured. For larger cloud column 
densities the disk may be gravitationally unstable and fragment into stars. Note that these 
limits on cloud column density are only a rough guide as the disk column density depends 
quadratically on the cloud speed and the angular momentum cancellation parameter (see eq 
7). 



4. Discussion 

AGN maser disks with accurate black hole mass determinations and Keplerian rotation 
curves sit well within the sphere of influence of their host black hole and are not self- 
gravitating. We showed that the outer radii of the megamaser disks in this sample scale 
linearly with host black hole mass. We then considered a scenario for the formation of these 
disks in which an incoming cloud temporarily engulfs the black hole and is partially captured 
(Sanders 1981; Bottema & Sanders 1986; Paper I). We showed that, for plausible estimates of 
the mass and angular momentum of the captured material, this process naturally reproduces 
the empirical linear relationship between maser disk size and black hole mass. The capture 
of clouds with column densities < 10^^'^ cm~^ results in a non-self-gravitating disks of the 
correct scale and sufficient column density to allow X-ray irradiation from the central source 
to reproduce thin megamaser disks. By contrast, the capture of a cloud with higher column 
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Fig. 2. — Constraints on the initial hydrogen column density A^h of the partially captured 
cloud as a function of the mass of the central black hole. Viable values he in the unshaded 
region - outside of this, the disk formed by capture is either unstable, is too massive, or has 
insufficient surface density for X-ray irradiation to produce the correct masing conditions 
(see text) . 
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density would instead create a self-gravitating disk giving rise to rapid star formation. In 

Paper I wc showed that this can explain the recent formation of the compact disks of stars 
within a fraction of a parsec of the Galactic center black hole Sgr A* (Paper I; Bonnell & 
Rice 2008; Mapelli et al. 2008; Alig et al. 2011). Note that the transition column density 
is sensitive to the speed of the incoming cloud and the degree of the angular momentum 
cancellation, though we have chosen plausible values of these parameters. 

This picture relies on the presence of dense clouds close to the nucleus of galaxies. 
The circumnuclear molecular ring (e.g. Christopher et al. 2005), 1.7pc from Sgr A* in 
our own Galaxy, and the circumnuclear rings found on scales of several parsecs from the 
center of numerous Seyfert galaxies suggest an ample supply of material. Recent simulations 
suggest that gas supply to galactic centers is controlled by angular momentum transfer 
from one massive gas clump to another during gravitational encounters (Namekata & Habe 
2011). Some of these inward-moving clouds may interact with their host supermassive black 
holes. The rate of migration of molecular material is estimated to give a potential black 
hole interaction rate of ~ 10^^ yr^^. Relating these estimates to the occurrence rate of 
megamaser disks requires estimates of the disk lifetimes and the near edge-on viewing angle 
needed for maser amplification. The single-dish detection rate of H2O megamasers in Seyfert 
2 galaxies and LINERs is about 15%, and about 20% of those have kinematics consistent 
with orbital motion on sub-parsec scales (Lo 2005). The estimated maser beaming angle, 
~ 10° (e.g. Maloney 2002) then imphes that ~ 20% of all Seyfert 2 galaxies possess similar 
disks. The VLBI follow-up results suggest that roughly half may be thin, Keplerian disks. 
The ordered kinematics of these disks suggest hfetimes of hundreds of orbital periods or 
more; the orbital period at 0.3 pc from a 10^ M© black hole is ~ 3 x 10^ yr implies lifetimes 
in excess of 10^ yr, consistent with the time scale needed to warp the disk via resonant 
relaxation (Alexander & Bregman 2009, 2011). The formation of gravitationally-unstable 
disks is likely to be as common simply because molecular clouds in the inner regions of 
galaxies tend to have column densities > lO^'^cm"^. These transient disks may also host 
megamasers (Milosavljevic & Loeb 2004), so that only a fraction of megamaser AGN may 
have disks that are very close to Keplerian enabling accurate black hole mass determinations. 

Finally, we note that the partial capture of a cloud imparts an impulse to the black 
hole / disk system of approximately M^iskV, where v ~ 200kms~^ is the incident cloud 
velocity. For plausible parameters this gives recoil velocities of ~ 10-20 km s~^; while this 
is potentially detectable, the recoil will be rapidly damped by dynamical friction on the 
surrounding stars, on a timescale of the stellar crossing time across the sphere of influence, 
~ 10^ years. 
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